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With the physical Higgs mass the Standard Model symmetry restoration phase transition is a 
smooth cross-over. We study the thermodynamics of the cross-over using numerical lattice Monte 
Carlo simulations of an effective SU(2)xU(l) gauge -|- Higgs theory, signihcantly improving on pre¬ 
viously published results. We measure the Higgs field expectation value, thermodynamic quantities 
like pressure, energy density, speed of sound and heat capacity, and screening masses associated 
with the Higgs and Z fields. While the cross-over is smooth, it is very well dehned with a width of 
only ~ 5 GeV. We measure the cross-over temperature from the maximum of the susceptibility of 
the Higgs condensate, with the result Tc = 159.5 ± 1.5 GeV. Outside of the narrow cross-over region 
the perturbative results agree well with uon-perturbative oues. 


I. INTRODUCTION 

The LHC particle accelerator has spectacularly con¬ 
firmed the Standard Model of the particle physics: the 
current combined results from ATLAS and CMS ex¬ 
periments point towards a Higgs boson with a mass 
125.1 ± 0.3 GeV [l|, and no direct experimental evidence 
of beyond the Standard Model physics has been observed. 

If the electroweak scale physics is fully included in the 
Standard Model, we are in the position to have a com¬ 
plete description of the high-temperature electroweak 
phase transition between the low-temperature broken 
phase, where the Higgs field has a non-vanishing expecta¬ 
tion value, and the high-temperature symmetric phase.^ 
The overall nature of the transition was settled already 
in 1995-98 using lattice simulations [m , which indicate 
a first-order phase transition for Higgs masses < 72 GeV, 
and a smooth cross-over otherwise. Most of these simula¬ 
tions were done using an effective 3-dimensional theory of 
the full Standard Model, derived using perturbation the¬ 
ory. The physics of the transition is non-perturbative due 
to the infrared singularities arising at momentum scales 
<g^T d, d|. These modes are fully captured in effec¬ 
tive 3-dimensional theories [IMl, which provide an eco¬ 
nomical and accurate way to study the non-perturbative 
physics at the cross-over. 

A cross-over means that the early Universe evolved 
smoothly from symmetric to broken phase without de¬ 
viating significantly from the thermal equilibrium, ren¬ 
dering e.g. electroweak baryogenesis [mil ineffective. 
Nevertheless, the cross-over could still influence other 
processes going on at the same time. An important ex¬ 
ample is the rate of B-l-L violation, sphaleron rate, which 
switches off somewhat below the cross-over temperature 


^ Although there is no real symmetry breaking phase transition, 
we use the conventional labels “broken” and “symmetric” to refer 
to the low- and high-temperature phases, respectively. 


m m and may affect the baryon number generation 
in leptogenesis . Decoupling of the dark matter 

may also be affected by the details of the Standard Model 
equation of state in certain scenarios [H [ll . 

The Standard Model cross-over with the physical Higgs 
mass of 125 GeV was recenty studied on the lattice in 
ref. m . The cross-over temperature was reported to be 
Tc = 159.5 ± 1.5 GeV. However, the focus of this study 
was the sphaleron rate as a function of the temperature, 
and for quantities related to thermodynamics the results 
were limited: as an example, the Higgs field expectation 
value ( 0 ^ 1 ^) was determined only at one lattice spacing, 
and the effects of the hypercharge U(l) held were ne¬ 
glected. 

Recently, Laine and Meyer [2^ used up to 3-loop per- 
turbative computations and existing results [ztL l28l [ to 
derive a relation between the trace anomaly of the Stan¬ 
dard Model and the Higgs held expectation value in the 
3-dimensional effective theory. Using the numerical re¬ 
sults from ref. [2^, they calculated several thermody¬ 
namic quantities across the cross-over. The combina¬ 
tion of perturbative calculations and lattice simulations 
avoids the infrared problems which make purely pertur¬ 
bative analysis unreliable near the cross-over. The set¬ 
up is related to the computation of the QCD pressure 
using effective theory |29| : however, for the Standard 
Model electroweak sector the couplings are smaller and 
the method can be expected to have much better accu¬ 
racy. 

In this work our goal is to improve on the state of the 
art in lattice simulations of the Standard Model cross¬ 
over: we include the U(l) gauge field, and we use sev¬ 
eral lattice spacings and large volumes, enabling us to 
do reliable continuum extrapolation. The model we sim¬ 
ulate is an effective three-dimensional theory containing 
SU(2)xU(l) gauge fields and a Higgs scalar. This partic¬ 
ular model has been used previously to study the Stan¬ 
dard Model phase transition [s^, also with an external 
hypermagnetic field [3l|. However, these older studies 
used unphysical Higgs masses. 
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As discussed above, measuring the Higgs condensate is 
essential for obtaining non-perturbative contributions to 
thermodynamic quantities. Thus, we pay special atten¬ 
tion to precise measurement of the gauge invariant expec¬ 
tation value and its susceptibility. The maximum 

of the susceptibility allows us to determine the pseud- 
ocritical temperature accurately, at = 159.6 ± 0.1 ± 
1.5 GeV. The first error is due to the precision of the lat¬ 
tice simulations and the second is a conservative estimate 
of the uncertainty of the effective theory description [T^ . 
The pseudocritical temperature is completely consistent 
with the result in ref. |^ . 

We use the condensate and the susceptibility to obtain 
several thermodynamic quantities across the cross-over: 
energy density, pressure, heat capacity, speed of sound, 
and equation of state parameter. For most quantities 
the magnitude of the effects of the cross-over are small, 
only at a percent level, but nevertheless clearly visible. 
We also determine the Higgs and screening masses 
and the — Z mixing. The emerging picture is fully 
consistent with a smooth and regular cross-over. The 
cross-over region, where observables deviate significantly 
from low- or high-temperature behaviour, is remarkably 
narrow, between 157 and 162 GeV. 

This paper is organized as follows. In section |IT] we 
describe the effective theory and in section IHII its imple¬ 
mentation on the lattice. The Higgs condensate and its 
susceptibility are discussed in section IlYl fundamental 
thermodynamic observables in section |V] and the screen¬ 
ing masses in section IVTl We conclude in section IVTIl 


Here Ai and Bi are the 3-dimensional SU(2) and U(l) 
gauge helds, and g'^ the dimensionful SU(2) and U(l) 
couplings, and </> a complex doublet. The SU(2)xU(l) 
local gauge transformation is 


The parameters appearing in the effective theory are di¬ 
mensionful. If we take one of the parameters, say g^, to 
set the scale, the dynamics then depends on the three 
dimensionless parameters x, y and z, defined as 
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The four parameters gf, g'^, A 3 and rug are definite per- 
turbatively computable functions of the Standard Model 
parameters {as{Mw), Gp, Mmggs, Mw, Mz, Mtop), and 
the temperature T. These have been computed through 
a set of 1 - and 2 -loop matching relations [l^, and are 
shown in figure [T] as functions of the temperature. The 
accuracy of the effective theory can be estimated to be 
at ^ 1% level, as discussed in refs. [H,[13- 

From figure [T] we see that only y has large temperature 
dependence. Indeed, from the effective theory point of 
view it is natural to choose y as the temperature variable, 
although we present our results in terms of the physical 
temperature. The transition is expected to happen near 
y = 0, which occurs at T = 162.1 GeV. 


III. LATTICE ACTION AND SIMULATIONS 


II. EFFECTIVE THREE-DIMENSIONAL 
DESCRIPTION 


At temperatures of order 100 GeV the gauge couplings 
in the Standard Model are small, and the Euclidean path 
integral contains a parametric hierarchy of energy scales: 
ttT, gT and g^T. While the harder scales can be reliably 
treated with perturbation theory, at fc ~ g^T we have 
to face the non-perturbative infrared physics Q. The 
non-perturbative physics can be captured into an effec¬ 
tive theory for the soft k ~ g^T scales, obtained by inte¬ 
grating over the harder scales using well-defined pertur¬ 
bative methods. This effective theory is purely bosonic 
and three-dimensional. The detailed description of the 
derivation of the theory can be found in refs. [ilH. 

The Lagrangian of the effective theory is the 3- 
dimensional SU(2)xU(l) gauge theory with a Higgs held 


L =-F“F“ 

4 y y 






where 


Fij = diAj - djAi - g^lAi, Aj], Ai = ^craA°; 

Bij = diBj - djBi 
Di = ^^ + igsAi + ig^Bil2. 


For the lattice implementation it is convenient to in¬ 
troduce a matrix parametrisation of the Higgs held by 
writing 



]_ { (t>2 ^ 

gl -4>l h J ■ 


( 5 ) 


Under an SU(2)xU(l) gauge transformation <i> trans¬ 
forms according to 

$(a:) ^ (6) 


The lattice action corresponding to the continuum theory 
0 is 

^ E E[i - + ^EE 

X i<.j X i<ij 

-PhEE (7) 

X i 

+ /32E5Tr$^(x)$(x) +/?4E 


Here the SU(2) and the (non-compact) U(l) plaquettes 
are 


( 2 ) 


P^jix) = Ui{x)Uj{x + i)u} {x + j)UJ {x), (8) 

aij{x) = ai{x) + aj{x + i) - ai{x + j) - aj{x). (9) 
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FIG. 1: The parameters of the effective theory 0 as functions 
of the physical temperature. 


The parameters of the lattice action (3c, I3 h, (^2 and Pa 
can be expressed in terms of the continuum parameters 
and lattice spacing a using relations [3l|, [s^ 


continuum quantity {((P p) in the MS scheme by 

(13) 

This relation again has corrections at order 0{a). The 
3-dimensional condensate is related to the physical Stan¬ 
dard Model Higgs condensate v by 

= 2{(jP(j))/T. (14) 

We use 3 lattice spacings, agl = 4,!Pc, with Pc = Q, 9 
and 16, with volumes up to (120a)^, as listed in table ID 
The measurements are done at 24 temperature values 
between 140 and 170 GeV. The simulation algorithm is 
a mixture of overrelaxation and heat bath update steps, 
described in detail (without a U(l) field) in [s^. The 
measurements of the observables are done after two rep¬ 
etitions of 4 full overrelaxation update sweeps followed 
by one heat bath update. The number of measurements 
varies between 10 000 and 350000, with most of the mea¬ 
surements done near the cross-over temperature. 


Pg = 4/5130 

volumes/a® 

6 

243 , 48® 

9 

24®, 32®, 60® 

16 

56®, 120® 


TABLE I: The lattice spacings and sizes used in the analy¬ 
sis. For each lattice, we use 24 temperature values between 
140 GeV and 170 GeV. The number of measurements varies 
between 10 000 - 400 000, with the largest number of mea¬ 
surements near the cross-over temperature. 
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where S = 3.17591. These relations become exact in 
the limit a —>■ 0 (Pc —t 00 ), with an error proportional 
to a. In refs. [13, [35l| partial 0(a) improvement of the 
relations was presented; however, for simplicity we do 
not implement this here. 

Finally, the gauge-invariant lattice observable 
(iTr4>l4>) is related to the renormalized 3-dimensional 


IV. HIGGS CONDENSATE AND THE 
PSEUDOCRITICAL TEMPERATURE 

The measurements of the Higgs condensate are 

extrapolated to the continuum with linear in a extrap¬ 
olation, as shown in figure [2] at a few selected temper¬ 
atures. For this quantity the finite volume effects are 
negligible, indistinguishable from the statistical accuracy 
of the measurements, and we use the largest volumes in 
our analysis. Below we analyze the finite volume effects 
of the susceptibility of p'^p in more detail and show that 
these also remain very small. Given the presence of the 
massless U(l) field, the smallness of the finite volume ef¬ 
fects is not a priori obvious, but is in agreement with the 
earlier results in [s^. 

The continuum extrapolation of (p^p) is shown in fig¬ 
ure [31 We compare the numerical result with the pertur¬ 
bative broken phase 2-loop Coleman-Weinberg commuta¬ 
tion [ 3 ^ and the symmetric phase 3-loop result [ 2 ^. In 
both cases the higher order contributions are estimated 
with shaded bands. The agreement between the pertur¬ 
bative results and the lattice results is remarkably good, 
especially in the symmetric phase where the perturbative 
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FIG. 2: The continuum limit of at a few selected tem¬ 

perature values. The statistical errors are too small to be 
visible at this scale. 



FIG. 3: The continuum result of compared with the 

perturbative broken and symmetric phase results. The shaded 
bands are estimations of unknown higher order corrections to 
perturbative results. The solid continuous line is an interpo¬ 
lation to the data. 


expansion converges quickly.^ There is only a narrow 
window of a few GeV around the cross-over temperature 
(corresponding to y ~ 0) where the perturbative expan¬ 
sions do not converge. 

The apparent good convergence in the symmetric 


^ Figure |3] can be compared with figure 2 in ref. [S^l . where the 
agreement between the lattice and the perturbative results is 
much weaker, due to the missing continuum limit of the lattice 
results. 



T/GeV 

FIG. 4: Above: susceptibility x<pf4, shown at /3 g = 6 , 9 and 
16, together with the interpolating functions. The continuum 
limit is shown with a heavy line. Below: As above, zoomed-in 
to the shaded band near the cross-over region. 


phase may be surprising, because in this phase the non- 
abelian gauge bosons are perturbatively massless, mak¬ 
ing the physics at soft momentum scales k ~ g^T non- 
perturbative Q. The excellent match between the lattice 
and the perturbation theory means that for the Higgs 
condensate their effect remains small. This can be con¬ 
trasted with e.g. the sphaleron rate, which is in essence 
completely determined by the soft physics. 

We define the pseudocritical temperature by the max¬ 
imum location of the dimensionless susceptibility 

= VT{[{^U)v - , (15) 

where {(j>'^(j>)v = 1/H J dVcj)^cj) is the volume average of 
This is shown in figure 01 for the largest simulation 
volumes at each lattice spacing. The use of the largest 
volumes is justified below. There is a well-defined peak 
near the cross-over temperature, however, the location 
of the peak has a clear lattice spacing dependence. Be¬ 
cause of the narrowness of the peak, the continuum limit 
extrapolation becomes delicate: at a fixed temperature, 
the values of at different lattice spacings have large 
and non-uniform variation, which can be clearly seen in 
the zoomed-in subplot in figure 01 Now a linear or a lin¬ 
ear + quadratic in a continuum extrapolation at fixed 
temperature does not give a reasonable result using the 
available lattice spacings. 

We obtain a much better controlled continuum limit if 
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FIG. 5: Continuum extrapolation of the maximum location 

of 

we first “undo” the shift in the temperature so that the 
peaks become aligned. This is perfectly compatible with 
0(a) lattice spacing effects, if the shift vanishes as 0(a) 
or faster as the continuum limit is approached. However, 
shifting the temperature requires that we interpolate the 
measurements of to continuous functions of T. We 
achieve this by fitting the measurements at each j3g to a 
differentiable interpolating function, XRt (T, /3 g)- The fit 
procedure and the choice of the function is described in 
appendix]^ The fitted functions at each /3g are included 
in figure 01 We note that commonly used reweighting 
methods do not work well in this case, because there is 
limited overlap in the probability distributions between 
neighbouring points over the large temperature range. 

We now define the pseudocritical temperature at each 
/3g by the location of the maxima of the fit functions 
/3g)- These are shown in figure O against 1/^g = 
ag§/4. The measurements fall on a remarkably straight 
line, allowing us to obtain the continuum limit pseudo¬ 
critical temperature, 

T’x»..(/3g)=T, + £, (16) 

with the result = 159.58 ± 0.06 GeV. The error bars 
conservatively include the results from both linear in a 
and linear -|- quadratic extrapolation. However, these er¬ 
rors are completely overwhelmed by the per cent-level 
uncertainties associated with the accuracy of the 3- 
dimensional effective theory description of the Standard 
Model. 

We use the above fit to shift the temperature variable 
of the functions Xfit as T — >■ T—CjP g- Now we obtain the 
continuum limit of x^t^ by making a linear extrapolation 
independently at each shifted value of T. The final result 
of the extrapolation is shown in figure 0] as a shaded band. 
The errors are propagated using the jackknife method. 
As can be observed, this method of extrapolation works 
very well. As an aside, if the fully 0(a) improved lattice- 
continuum relations corresponding to equations (iinHni) 
were available, the 0{a) shift in T should be cancelled 



T/GeV 

FIG. 6: The finite volume behaviour of the susceptibility 
X^i^j shown for the lattices of size Lja = 24, 32 and 60 at 
Pa = 9. Only the smallest volume shows appreciable devia¬ 
tion, which is much below the lattice spacing effects in figure 

a 


automatically. 

Finally, in figure [6] we show the behaviour of the sus¬ 
ceptibility at different volumes at /3 g = 9. In this case 
we have 3 volumes available, with physical linear lattice 
sizes L Ri 27/T, 36/r and 68/T (at Lja = 24, 32 and 
60, respectively). The larger volumes agree within sta¬ 
tistical errors, and the smallest volume deviates slightly 
only at the susceptibility peak. Compatible behaviour 
is seen at other lattice spacings, where the smallest vol¬ 
umes have L>36/r. Thus, we conclude that the largest 
volumes have negligible finite volume corrections for the 
Higgs condensate and its susceptibility. 

In the absence of massless modes the finite volume ef¬ 
fects vanish exponentially in L. In this case the system 
has a massless U(I) gauge boson, and thus one could 
expect to observe power law finite volume effects. How¬ 
ever, the hypercharge U(l) field participates weakly to 
the transition, only by mixing with by the Weinberg 
angle. In part, this explains the smallness of the finite 
volume effects.^ Similar behaviour has been observed in 
earlier lattice studies 

In order to facilitate the detailed analysis of ther¬ 
modynamic functions we also interpolate the 
measurements. We do this by interpolating at each lat¬ 
tice spacing using a spline fit, shifting the temperatures 
by the same amount as for x<?it 0 and extrapolate to con¬ 
tinuum. The result is included in figure[3]as a continuous 
line. The temperature shift before the continuum extrap¬ 
olation changes the result by a very small amount which 
not distinguishable in the figure. For possible future use 
the interpolated results for p'^cj) and x<^t,^, together with 
other thermodynamic measurements, can be downloaded 
from [ 3 ^. 


® Strong external hy per magnetic fields do nevertheless affect the 
transition, see e.g. l3ll . 
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FIG. 7: The interaction measure A = (e — 3p)/T‘^ (top left); energy density e and pressure p (top right); heat capacity 
Cv = de/dT (bottom left); speed of sound squared and the equation of state parameter w — pje (bottom right). The error 
bands are a combination of the statistical errors and renormalizaton scale variation p. — (0.5 ... 2)irT. The energy density and 
the pressure are affected by a systematic uncertainty of order 1%, indicated with a vertical arrow. 


V. THERMODYNAMICS OF THE 
CROSS-OVER 


Recently, Laine and Meyer [2^ showed how one can 
combine perturbative calculations and effective theory 
simulation to obtain the Standard Model pressure, en¬ 
ergy density and other thermodynamic quantities derived 
from these. As an input they used the simulation results 
from ref. . Because these results use only a single lat¬ 
tice spacing and ignore the hypercharge U(l), we revisit 
the calculation using our improved data. 

The fundamental thermodynamic quantitity is the in¬ 
teraction measure (“trace anomaly”) 


^ _ e{T) - 3p{T) d p{T) 
Ti dT T4 ■ 


(17) 


In ref. the interaction measure is split into three 
parts: A = Ai-t-A 2 -l-A 3 , where Ai includes effects from 
breaking of scale invariance by quantum corrections, A 2 
effects from the Higgs condensate, and A 3 comes from 
vacuum subtraction. For our purposes it is convenient to 
express it as 


A(T,/ 2 )=A(r;/i) + H(T;/ 2 ) 


{4>^4>){T) 


T 


(18) 


where p is the MS renormalization scale and the functions 
A and B can be computed following ref. [ 2 ^ , giving A up 


to parametric order g^. However, here we use the 0{g'^) 
expression for the function A, because it has been argued 
that the 0{g^) contribution, which is mostly due to QCD 
contributions, leads to an underestimate of the pressure 
and the energy density [ 2 ^ . The resulting A is shown in 
figure [T] 

From figure [3] we can see that the direct perturba¬ 
tive computation of {(j)'(j)){T) fails to converge near the 
cross-over temperature. On the other hand, the func¬ 
tions A{T-,p,) and B(T;p,) do not suffer from this prob¬ 
lem. Therefore, by measuring the Higgs condensate non- 
perturbatively on the lattice, we obtain well-defined ex¬ 
pressions for thermodynamic quantities across the cross¬ 
over. 

The pressure is obtained from A by integration: 


r4 ^ T' 


(19) 


In order to evaluate this we need to fix the pressure 
at a reference temperature Tq. We use the results in 
refs. [ 2 ^ 113 at the lowest temperature in our temper¬ 
ature range, Tq = 140GeV: p{To)/Tq = 11.173. The 
estimated uncertainty in this value is of order 1%; we 
discuss this together with other systematic errors at the 
end of this section. 

From A(T) and p{T) we can obtain other thermody- 
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namic functions: the energy density e/T^ = A + 
the entropy density s = p' = {e+p)/T, the heat capacity 
Cv/T^ = e'/T^ = 7A + TA' + Up/T'^, the speed of the 
sound squared = p' je' and the equation of the state 
parameter w = pje. The quantities involving A' can in 
principle be calculated directly from (jl8|) by numerical 
differentiation. However, we obtain a numerically more 
robust result by realizing that to the parametric accuracy 
needed 


d _ dy d ^ dy 

dT T ^ dTdy T ~ 


The direct numerical differentiation of (j)) gives com¬ 
parable results, but because it was obtained through con¬ 
tinuum extraplation of spline interpolations the deriva¬ 
tive of it becomes somewhat jagged. 

The final results are collected in figure [3 with er¬ 
ror bands which combine quadratically the statistical 
errors and the variation in the renormalization scale 
p, = (0.5T... 2)ttT. The renormalization scale dominates 
the error bands except in quantities where A' contributes, 
where the statistical errors are of comparable magnitude. 
The results show similar features than the earlier ones 
published in ref. , but with improved accuracy and 
reliability, because the Higgs condensate and suscepti¬ 
bility have been properly extrapolated to the continuum. 
The tabulated numerical results can be downloaded from 



Figure [3 clearly shows that the cross-over region, de¬ 
fined as the temperature range where the thermodynamic 
variables deviate from their (perturbative) hot or cold 
phase behaviour, is narrow, between ~ 157 and 162 GeV. 
This is precisely the region where perturbation theory 
does not converge and lattice simulations are needed. 
The variation of the thermodynamic quantities across the 
cross-over is very mild, only at few per cent level. This 
is a natural consequence of the fact that only a couple of 
degrees of freedom of the Standard Model are strongly 
sensitive to the cross-over. The behaviour is markedly 
different from the QCD cross-over, see e.g. [s^ for a re¬ 
cent review. 

In addition to the statistical errors and renormaliza¬ 
tion scale ambiguity shown in figure [3 there are theoret¬ 
ical uncertainties that can be estimated to be at percent 
level [1^, especially affecting the energy density and the 
pressure. As an indication of these a 1% bar is shown 
in figure [3 and it clearly dominates over the error bands 
of e and p. This uncertainty mostly pertains to absolute 
normalization in these quantities, whereas the shape of 
the curves can be expected to be preserved to much bet¬ 
ter accuracy. In particular, the characteristic features of 
the thermodynamic functions around the cross-over are 
mostly due to the evolution of {(j)' (j)) and hence more ac¬ 
curately determined. 
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FIG. 8: Higgs and screening masses, measured from 60®, 
Pg = 9 and 120®, /3 g = 16 lattices. 


VI. SCREENING MASSES 

The screening masses of the Higgs and SU(2) gauge 
fields are direct probes of the “softening” of the physics 
at the cross-over. We measure the Higgs, IF® and hyper¬ 
charge B spatial correlation functions and extract the 
screening masses (inverse correlation lengths) as func¬ 
tions of the temperature. The Higgs and IF® correla¬ 
tion functions are measured with lattice operators corre¬ 
sponding to H = Tr<I>l$ and IF® = Tr pro¬ 

jected to zero momentum. For example, the Higgs cor¬ 
relation function is (here in units of lattice spacing a): 


^h{z) = ^ X! {H(xi,yi,z')H{x 2 ,y 2 ,z' + z)), (21) 

where Xi, t/i, z, z' are {x, y, z) coordinates compatible with 
the periodic boundary conditions. The local operators 
corresponding to IF®' and IF“ (linear combinations of 
Tr with a = 1,2) are not gauge invariant under 

U(l) and we do not measure them here. To enhance the 
signal, we use lattice operators which are smeared and 
blocked recursively 3 times, as described in ref. [s^. The 
screening masses extracted from the correlation functions 
are shown in figure [51 We use the largest lattices at 
/3 g = 9 and /3 g = 16, with results which are consistent 
within the statistical errors. 

In the broken phase, where has a large expec¬ 
tation value, the operators couple appropriately to the 
physical Higgs and Z-boson excitations. At small tem¬ 
peratures the Higgs screening mass is larger than that of 
IF®, whereas near the cross-over temperature the Higgs 




























becomes sensitive to the near-critical behaviour and the 
mass becomes small but still remains non-zero. Even at 
its largest, the Higgs correlation length is smaller than 
10/T, which is substantially smaller than the largest lat¬ 
tice sizes ^ 70-80/r. 

In the symmetric phase, the non-abelian gauge fields 
are confining, and the operators couple to bound states 
of two scalars. The correlation functions become noisy 
and the screening masses increase rapidly. 

The U(l) gauge field correlation function can be used 
to measure the ^-Z mixing, i.e. the effective Weinberg 
angle. We define the operator 

= X! CKi2(xi,X2,z)e‘‘P'^, (22) 

Xi ,X2 

where the sum is taken over the plane (xi,X 2 ), ctij is 
the (non-compact) hypercharge U(l) plaquette ([9]) and p 
is a transverse momentum vector compatible with peri¬ 
odic boundary conditions: {pi,P 2 ,P 3 ) = 2 Tr/N{ni,n 2 , 0 ) 
with integer ni and n 2 - In our measurements we use 
the smallest non-vanishing momentum, with |p| = 2'ir/N. 
At p = 0 the operator Op vanishes, due to the periodic 
boundary conditions. The correlation function 

= + (23) 

t 
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FIG. 9: The effective 7 —Z mixing as a function of the temper¬ 
ature. The dashed lines show the 1-loop perturbative results. 


VII. CONCLUSIONS 


has the long distance behaviour 


G{z) 




ap 


Wg 




(24) 


where is the photon screening mass and gives 
the projection of the operator to the hypercharge U(l) 
field, in effect yielding the temperature-dependent effec¬ 
tive mixing angle. At tree level, A~f = 1 in the symmetric 
phase and A^ = cos^ 9w in the broken phase. 

The photon screening mass vanishes within our 
measurement accuracy at all temperatures. The projec¬ 
tion Aj is shown in figure [5] for /3 g = 9, 60^ lattice. The 
measurement is noisy, but we can observe that A^ Ri 1 in 
the symmetric phase down to the cross-over temperature, 
and it starts to decrase as the Higgs field expectation 
value grows at lower temperatures, slowly approaching 
the tree-level value. 

Beyond tree-level perturbative estimates for the be¬ 
haviour of A^ can be obtained by calculating at 1-loop 
order the residue of the 1/fc^ pole in the {BiBj) corre¬ 
lator. In the symmetric and broken phases one obtains 


^symm. _ ^ _ 

^broken ^ 


11 gj sin^ 

12 TTTTH 


(25) 

(26) 


where mw is the perturbative W mass. These expres¬ 
sions clearly anticipate the behaviour we observe on the 
lattice, although they diverge as y —>■ 0±. 


We have accurately determined the Higgs field expec¬ 
tation value and its susceptibility across the Standard 
Model cross-over using lattice simulations of an effective 
3-dimensional theory. Defining the cross-over temper¬ 
ature by the maximum of the susceptibility, we obtain 
Tc = 159.6 ± 0.1 ± 1.5 GeV, where the first error is due 
to the statistical accuracy of the lattice computation and 
the second one is the estimated uncertainty of the effec¬ 
tive theory approach [I^ [2^ . Following the approach of 
Laine and Meyer [2^, these results were used to obtain 
the behaviour of basic thermodynamic quantities, includ¬ 
ing energy density, pressure, heat capacity and the speed 
of sound, across the cross-over. There is a well-defined 
cross-over region where thermodynamic quantities devi¬ 
ate from the low- or high-temperature behaviour. This 
region is quite narrow, between 157 and 162 GeV. The re¬ 
sults are consistent with the standard picture of the elec- 
troweak cross-over: Higgs and W modes become softer 
but not critical, and the U(l) field remains massless at 
all temperatures. 

Overall our results are compatible with the analysis in 
ref. [ 2 ^ using lattice data from ref. [l^. Howeever, our 
results are significantly improved numerically: we have 
much larger volumes with higher statistical accuracy, the 
data is extrapolated to the continuum and we include 
the U(l) field in the effective theory. Thus, our results 
form an important consistency and reliability check of 
the earlier results. 

For phenomenological applications the thermodynamic 
quantities here can be combined with existing low- 
and high-temperature perturbative results. This has 
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been done in ref. and we do not repeat this analysis 
here. 

The errors are dominated by systematic uncertainties 
at per cent level, but since these are expected to be 
roughly constant at the temperature range considered, 
they would mostly shift the scale of the results without 
changing their behaviour. The cross-over is clearly visi¬ 
ble in distinct features in the heat capacity, the speed of 
sound and the interaction measure. The relative magni¬ 
tude of these features is small, only at few % level, which 
already follows from the fact that only few degrees of 
freedom of the Standard Model are sensitive to the early 
stages of the transition. 

In ref. [ 2 ^ the central value of the cross-over tempera¬ 
ture was reported to be 159.5 GeV, in practice equal to 
our continuum result here. This work was done without 
the U(l) held and only at one lattice spacing, /3 g = 9, 
using relatively small statistics. However, one cannot 
make the conclusion that the effects of the U(l) held 
can be neglected in the transition; from hgure [5] we see 
that at / 3(3 = 9 the maximum of the susceptibility is at 
~ 159.0 GeV, which is a numerically signihcant deviation 
from the continuum result. Thus, the apparent excellent 
agreement is partly coincidental, due to a different dehni- 
tion of the cross-over temperature and limited accuracy 
in [ 2 ^. Nevertheless, we can estimate that the presence 
of the U(l) held in the effective theory shifts the cross¬ 
over temperature <0.5 GeV. 
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Appendix A: Interpolating function for the 
susceptibility 

In this subsection we describe the method used for the 
interpolation of the susceptibility measurement at differ¬ 


ent lattice volumes and /Sc-values in section El Instead 
of the temperature, it is convenient to do the interpola¬ 
tion in terms of the 3-dimensional variable y = Ta\lg^, 
which is actually the variable used internally throughout 
our analysis. 

The functional form of the ht is motivated by the 
sharp susceptibility peak near y = 0 {T = 162.1 GeV), 
approaching almost flat behaviour away from the peak. 
Thus, we construct a differentiable function which ap¬ 
proaches power series in 1/y on both sides of the cross¬ 
over. For our data we obtain very good results with the 
function 

=- E^iy) + E,iy) -’ 

(Al) 

where 

ys=y-ai 

s = 1 /V {yl + 0,2 

E^{y) = l/E2{y) = e-^^y^ 

Here oi... 03 , cq ... and do ... dn 2 are ht parameters, 
s corresponds to a regulated |?/|“^ and the exponential 
factors Ei ^2 are introduced to smoothly switch between 
the c-series on the low temperature side and the d-series 
on the high temperature side of the cross-over. Param¬ 
eter oi corresponds roughly to the location of the peak, 
and 02 > 0 regulates the singular behaviour at ys 0 . 

Using ni = 3, 712 = 2 we obtain a good ht to all mea¬ 
surements (x^/d.o.f. = 13/16, 20/16 and 10/16 for the 
largest volumes at /3 g = 6 , 9 and 16, respectively), as 
shown in hgure [H Adding higher order terms does not 
improve the ht, and the hnal function remains practi¬ 
cally identical, with variation remaining well within the 
statistical errors. 

The ht function is not unique, and the ht parame¬ 
ters have no direct physical meaning. The function only 
serves to interpolate the measurements and cannot be ex¬ 
tended beyond the range of data shown here. Naturally, 
signihcantly more accurate data or extended temperature 
range would require a more general ht function or some 
other means of interpolation. The error propagation is 
done with jackknife analysis. 
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